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The Fano factor, an observable quantifying fluctuations of product generation by a single enzyme, 
can reveal information about the underlying reaction scheme. A lower bound on this Fano factor 
that depends on the thermodynamic affinity driving the transformation from substrate to product 
constrains the number of intermediate states of an enzymatic cycle. So far, this bound has been 
proven only for a unicyclic network of states. We show that the bound can be extended to arbitrary 
multicyclic networks, with the Fano factor constraining the largest value of the effective length, 
which is the ratio between the number of states and the number of products, among all cycles. 


I. INTRODUCTION 

The Michaelis-Menten (MM) expression for the rate of 
product formation in terms of the substrate concentra¬ 
tion is a central paradigm in enzymatic kinetics [I|- In 
this case, the reaction scheme is very simple, with one 
intermediate state, corresponding to a compound of en¬ 
zyme and substrate, and the further assumption that the 
generation of product is irreversible. 

However, enzymatic reactions schemes are generally 
more complex. They can involve several intermediate 
states and, possibly, different enzymatic cycles, arising, 
for example, from different conformational forms of the 
enzyme. Single molecule experiments yield precise 
data on the dynamics of single enzymes, including fluctu¬ 
ations related to the rate of product formation. Obtain¬ 
ing information about the underlying enzymatic network 
through the measurement of fluctuations constitutes a 
promising field known as statistical kinetics [M3- 

One important quantity characterizing the dispersion 
of a probability distribution is the Fano factor [l0l - ll3j . 
which is the variance squared divided by the mean. For 
unicyclic enzymatic networks the Fano factor related to 
the fluctuating product creation is bounded from be¬ 
low hy the inverse of the number of intermediate states 
[ICi [ 14 L [ 1 ^ . Hence, its measurement provides informa¬ 
tion on the minimal number of intermediate states in the 
enzymatic cycle. 

This bound on the Fano factor does not take into ac¬ 
count the affinity driving the enzymatic process, which is 
given by the chemical potential difference between sub¬ 
strate and product. Often the chemical reaction leading 
to product generation is assumed to be irreversible, as 
is the case of the MM scheme, which would lead to a 
formally divergent affinity. Considering only thermody¬ 
namically consistent models of enzymatic kinetics, with 
no irreversible transitions, we have recently shown that 
a stronger bound on the Fano factor that depends on the 
affinity can be obtained [l^ . 

Proofs of these bounds assume a unicyclic network of 
states [3;[I3- However, chemical reaction networks can 
involve many different cycles. In this paper we obtain an 
affinity dependent bound for arbitrary multicyclic net¬ 
works of states. Calculating the Fano factor for different 


multicyclic networks exactly we observe that a version 
of the bound for the unicyclic case can be adapted to 
general networks. 

The paper is organized as follows. In section |H1 we 
state the bound on the Fano factor, which can be proved 
for unicyclic networks. In section nil we infer from ana¬ 
lyzing two case studies that this bound can be adapted to 
multicyclic networks. The constraints imposed by mea¬ 
surements of the Fano factor on network topology for the 
case where the affinity is known are discussed in section 
m We conclude in section lYl The method used to cal¬ 
culate the Fano factor is explained in Appendix In 
Appendix [B] we provide explict calculations for a sim¬ 
ple example. Appendix [C] contains two further examples 
of multicyclic networks demonstrating the validity of the 
bound. 


II. UNICYCLIC NETWORKS 

A. MM kinetics 

A standard example of a unicyclic network of states is 
provided by MM kinetics with reversible product gener¬ 
ation, which is represented as 

k'^ 

E + S^ES^E + P, ( 1 ) 

k- k- 

where E is the enzyme, S a substrate molecule and P a 
product molecule [l|. The transition rate is propor¬ 
tional to substrate concentration [S'], whereas the tran¬ 
sition rate is proportional to product concentration 
[Pj. It is assumed that [S] and [P] are large enough so 
that they remain approximately constant. Moreover, the 
difference in chemical potentials drives the for¬ 

mation of product creation with an affinity 

jA = fig fip. ( 2 ) 

In terms of the transition rates, the affinity is given by 
the generalized detailed balance condition [ 13 , [T^j 

A/iksT) = ln[ktk+/ik^k^)], (3) 


2 


where T is the temperature and ks Boltzmann’s con¬ 
stant. In the following, we set ksT = 1, so that the 
affinity and chemical potentials are measured in units of 
keT. 

For a single enzyme, the time evolution of the model 
represented in Eq. m is determined by the master equa¬ 
tion, which reads 


= (fc+ + k-)P2{t) - {k+ + k^)P,{t), (4) 


where Pi(t) is the probability that the enzyme is free, 
i.e. in state E, at time t and P 2 {t) = 1 — Pi{t) is the 
probability of state ES, with the enzyme bound by a 
substrate. Independent of initial conditions, the sta¬ 
tionary probability is reached for large t and given by 
Pi = -I- k^)/(k^ Pk'l -I- -I- fc^), which is obtained 

by setting = 0. In this stationary state, the rate at 
which the product P is generated is given by 


V = (1 - Pi)k^ - Pik2 


^•+ 4 -+ — h~h~ 

k'^ -\- k^ -\- k-^ F ^^2 


(5) 



FIG. 1: Enzymatic network for the four-state model. The 
links represented by a are related to chemical reactions in¬ 
volving substrate S and the links represented by b are related 
to chemical reactions involving product P, in agreement with 
Eqs. (1131) and (11411 . The dashed links form a four-state cycle 
that leads to the minimal Fano factor F. 


then we would get TV = 3. Furthermore, n is the number 
of consumed substrate molecules, which is equal to the 
number of generated product molecules, in the cycle. The 
affinity of the cycle then becomes 

n^ = ln(r+/r_), (11) 


We are interested in a random variable X that counts 
the number of generated product molecules. Moreover, 
we consider that t is large enough so that 


{X) = vt, 

(6) 

where the brackets denote an average over 
realizations. The Fano factor is defined as 

fluctuating 

{X^)-{Xf 

(7) 

{X) F’ 

where 


D = ((X2) - {Xf)/(2t) 

(8) 


is known as dispersion. Whereas the rate of product cre¬ 
ation V has a simple expression in terms of the stationary 
probability, this is, in general, not the case for the dis¬ 
persion D 113. For the model in Eq. o the dispersion 
becomes [Ill3 


D = 


ki ^2 +^ 1^2 “ 


2{k 


k+ 


fci 




(9) 


where V is given in Eq. ®. 


where r+ = k:^, F. = k~, and A is given by 

Eq. ([2]). The Fano factor for this general case is defined 
as in Eq. Q. 

An important bound valid for unicyclic machines is 
P>n/N [13,Hd. Hence, measurements of the Fano 
factor provide a bound on the number of intermediate 
states N. Recently, we have obtained a stronger, affinity 
dependent bound on F given by (l3 | 

This improved bound becomes F > n/N if A is large. 
An important point for the discussion in the next section 
is that this bound is saturated if the rates are uniform, 
i.e., kf = F^^'^ and k~ = for i = 1 ,..., Af (T3| - 


III. MULTICYCLIC NETWORKS 


Enzymatic networks can be more complex than a single 
cycle. In this section we address the question of how 
the bound m can be adapted to arbitrary multicyclic 
networks. 


B. Bound on F for unicyclic networks 

We consider a general unicyclic network with N states 
for the enzyme. For example, for the model in Eq. o 
N = 2. If we had included the intermediate state EP, 
leading to the reaction scheme 

fc'*” k'^ 

E + S ^ ES ^ EP E + P, (10) 

^2 ^3 


A. First case study 

As a first example of a multicyclic network of states 
we consider an enzyme E that can also be in a second 
conformational state E*, see Fig. [T] The effect of having 
such conformational states on the MM expression has 
been recently analyzed in [lllll^. There are four possible 
states, which are identified as E = 1, ES = 2, E* = 3 


















3 


and E*S = 4. Besides the reaction scheme in Eq. o, 
which now takes place with rates 

k-t n kn-t 

E + S^ES^E + P, (13) 

^21 ^^2 

there are also the reactions 

E* + S ^ E*S ^E* + P, (14) 

^43 fc34 



^31 


and 



ki2 


E*S. 


(15) 


(16) 


The enzymatic network corresponding to these chemical 
reactions is shown in Fig. [TJ 

The transition rates for this model are not all inde¬ 
pendent. For thermodynamic consistency, the following 
constraints must be fulfilled. Since the affinity of cycles 
where one S is consumed and one P generated is A, from 
the cycle in Eq. m we get 


In 


ya yb 

_ A 

ya yb ~ 


and, from the cycle in Eq. (HI we obtain 


(17) 


In 


ua ub 
^34^43 

ha hb 
^43^34 


= A. 


(18) 


Moreover, the cycle 


E + S 


k-i o 


ES ^ E*S ^ E* 


koi 


■S^E + S, (19) 

ki3 


which does not lead to product formation, has affinity 0, 
leading to the constraint 


In ^12^24^43^31 


( 20 ) 


Hence, these 3 constraints for 12 rates lead to 9 indepen¬ 
dent transition rates. 

Denoting the probability of being in state i at time t 
by Pi(t), the master equation reads 


jPit) = LPit) ( 21 ) 

where P(t) = {Pi{t), P2{t), Pait), P4{t)} and the stochas¬ 
tic matrix L is 


/ 

-ri 

ya 1 yb 

"21 ' "'21 

ksi 

0 \ 


ya 1 yb 

"'12 w n/42 

-r2 

0 

k42 


kl3 

0 

-rs 

ha 1 hb 
^43 ^ /^43 

V 

0 

k24 

ya 1 yb 
"-34 w /i34 

-r4 / 



A 


FIG. 2: The Fano factor E as a function of the chemical po¬ 
tential difference M, evaluated at randomly chosen rates, for 
the model shown in Figure [T] The red solid line corresponds 
to the bound (1/4) coth(yl/ 8 ). The transition rates fc 2 i, ^ 43 , 
and fci 3 are eliminated with the relations ^21 = e'^fci 2 fc 2 i/fci 2 , 

k4$ = G'^k^4k4;^lk^4, and kl3 = ki2k24k43k3-i/(^k2lk42k34'). 
The remaining transition rates are chosen by generating a 
random nnmber y between —5 and 5 and then taking k = 10^. 
This symmetric choice of the rates lead to the property that 
only a few points go below (1/2) coth(^/4), which is shown 
as a blue dashed line, indicating that the region in the space 
of transition rates where F is close to the bound is small. 

where ri = fc “2 + ^12 + ^i3, = ^21 + ^21 + ^ 24 , 

= ksi + k ^4 + k^ 4 , and r4 = fcfg -|- k^^ + k42. 
The stationary probability vector P is obtained from 
LP = 0. Furthermore, the rate of product generation 
is E = P 2 ^ 2 i ~ Piki 2 + ^4^43 — ^3^34- The calculation of 
the dispersion D, as explained in Appendix 1X1 leads to a 
quite long expression for the Fano factor E as a function 
of the affinity A and 9 independent transition rates. 

In Fig. O where we evaluate the Fano factor at ran¬ 
domly chosen transition rates obeying the constraints, we 
see that F is bounded by the right hand side of Eq. m 
with N = A and n = 1. This bound is related to cycles 
in the network of Fig. [2] which have N = A states. An 
example of such a cycle, as represented by the dashed 
links in Fig. [TJ is 

E* + S^E*S^ES$^E + P^E* + P (23) 

fcJs fc24 k^^ ^31 

Choosing uniform transition rates for this cycle and con¬ 
sidering the transition rates that are not part of the cy¬ 
cle to be small, F should indeed reach the bound in Eq. 
(IT^ with N = A. For example, setting ^34 = ^42 = 
^21 = ^13 = exp(A/4), fcfg = k24 = k^^ = ksi = 1, 
^12 = ^43 = ^exp(3M/4), and k^i = k^^ = k in the 
expression for F, by taking the limit fc —^ 0 we obtain 
analytically 

i^=icoth(:^). (24) 

In Appendix jBj we consider a simpler network, where a 
similar calculation can be performed explicitly. As an 
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E=1 



FIG. 3: Enzymatic network for the six-state model. The 
dashed links are part of a five-state cycle leading to the min¬ 
imal Fano factor. 

independent check that the Fano factor is indeed larger 
than the above expression we have also minimized it nu¬ 
merically. Hence, even though the MM reaction scheme 
in Eqs. m and (fni) has two states, the Fano factor for 
the full model in Fig. [I]can be smaller than the bound in 
Eq. (fT^ with N = 2: the conformational changes in the 
enzyme allow for larger cycles with = 4. The bound 
in Eq. m for unicyclic networks thus also applies for 
the multicyclic network of Fig. [T] with N being the size 
of the largest cycle in the network. 



A 


FIG. 4: The Fano factor F as a function of the chemical 
potential difference A, evaluated at randomly chosen rates, 
for the model of Fig. [S] The red solid line corresponds to the 
bound (1/5) coth(^/10). From the constraints in Eq. (1291) 
we eliminate the transition rates fcai = e^fe 2 i 1532 ^ 13 /(^ 12 ^ 23 ), 
fees = s'^kieksske^/^ksskss), 1552 = e'^fess 15251532 /( 15351523 ), and 
I 545 = I542l554l535l523/(l524l553l532). The remaining transition 
rates are chosen by generating a random number y between 
—5 and 5 and then taking k — 10^. The fact that the points 
are far from the bound indicates that the region in the space 
of transition rates where F becomes close to the bound is 
small. 


B. Second case study 


As a second example of a multicyclic network we con¬ 
sider an enzyme that can bind two substrate molecules 
[l|, as shown in Fig. [3] Intermediate states with an 
enzyme bound with product molecules are also allowed. 
The model has 6 states which are identified as E = 1, 
ES = 2, EP = 3, ESS = 4, ESP = 5, and EPP = 6 . 
The possible reactions are 

E + S ^ ES ^ EP ^ E + P, (25) 

^21 ^32 ki3 


s ^ ESS ^ ESP ^ES + P, (26) 

A;42 k54 k25 


lpqu 


(27) 


There is also another three-state cycle, which is the mid¬ 
dle cycle in Fig. [3] that does not involve further reactions. 
The stochastic matrix L, defined in Eq. (EU, for this 
model is 


/ 

-ri 

^21 

kei 

0 

0 



ki2 

-r2 

^32 

^42 

ke2 

0 


ki3 

^23 

-re 

0 

^53 

kes 


0 

k2i 

0 

- 4-4 

keA 

0 


0 

^25 

^35 

^45 

-re 

kee 

V 

0 

0 

^36 

0 

kee 

-re / 


(28) 


where ri = ki 2 + ki3, r 2 = fc 2 i+ 1523 + ^24 + 1525 , rs = I 531 + 
I 532 + I 535 + 1536, ?’4 = I 542 + I 545 , r5 = I 552 + I 553 + I 554 + 1556, 
and re = kes + fees- The 18 transition rates fulfill the 4 
constraints 

A _ In ^12^23^31 _ ]y. k 24 k 4 ?,kr ,2 

^21^32^13 ^42^54^25 

_ In ^35^56^63 _ k23k35k52 /'9Q'l 

^53^65^136 ^32^53^25 ’ ^ ^ 

which come from the fact that the affinity for each of the 
four three-state cycles indicated in Fig. |3]is A. The rate 
at which product is generated can be written as E = 
Psksi — Pikis + Pekes ~ Pakse + Pbke 2 — P2k25- The 
calculation of the dispersion D, which is explained in 
Appendix leads to an even longer expression for P as 
a function of A and 14 independent transition rates. 

The largest cycle for the network in Fig. [3] is the 
six-state external cycle going through all states. Even 
though N = 6 for this cycle, the number of consumed 
substrate and generated product is n = 2. Hence, the 
bound on the Fano factor in Eq. dm for this cycle with 
N = 6 and n = 2 is the same as the bound for the 
four cycles with N = 3 and n = 1 indicated in Fig. |31 
There are also cycles with N = 5 states and n = 1, 
with an example indicated by the dashed links in Fig. 
|3l The latter cycle gives the bound on the Fano factor, 
which is reached when the multicyclic network becomes 
effectively this cycle with uniform rates. This limit can 
be achieved by setting the transition rates, e.g., in the 
following way. For the five-state cycle we take the uni¬ 
form rates ^32 = ^24 = ^45 = kee = kee = exp(A/5) 
and ^23 = ki 2 = ke 4 = kee = fcse = 1- As the tran¬ 
sitions related to the links EP — ESP and ES — ESP 
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should be slow, we set ^35 = k ^2 = fcexp(3.4./5) and 
^53 = ^25 = k, with k being small. The transition to 
leave the cycle to state E should also be slow, whereas 
the transitions to return to the cycle from state E should 
be fast, as with ki 2 = exp(3yl/5)/fc, ^21 = k, ki^ = 1/k, 
and ^31 = A:exp(3^/5). With this choice for the tran¬ 
sition rates in the expression for E, by taking the limit 
k ^ 0, we obtain analytically 

We verified with numerical minimization and evaluation 
at randomly chosen transition rates, as illustrated in Fig. 
m that the Fano factor does not go below this bound. 


C. Conjecture for arbitrary networks 

The results obtained with the two case studies and the 
two further examples shown in Appendix [C] lead us to the 
following conjecture. Given an enzymatic process where 
a substrate S is consumed and product P is produced 
with an arbitrary enzymatic network, the Fano factor is 
bounded from below by 

F > — coth ( ^ , (31) 

- Mmax V2M„iax/ 

where M^ax is the largest value of the effective length 
M = N/n, i.e., the length per product molecule, among 
all cycles in the multicyclic network. A posteriori, the 
idea behind the conjecture is simple: a multicyclic net¬ 
work cannot have a Fano factor smaller than the one 
of a unicyclic network with the largest effective length 
A/max- If the rates are chosen such that the network is 
fully dominated by this cycle with uniform rates, chang¬ 
ing the rates will add cycles with smaller M leading to 
an increase in F. 

A somewhat related result that for arbitrary Markov 
processes the relative uncertainty associated with the 
probability of the time to reach an absorbing state is 
bounded by the inverse of the total number of states in 
the network has been proven rigorously [ 2 ^. This rela¬ 
tive uncertainty is equal to the Fano factor for the special 
case of a unicyclic network with at least one irreversible 
transition | 10 l| . which amounts to a formally divergent 
affinity. 


IV. FANO FACTOR AS A DIAGNOSTIC TOOL 
FOR NETWORK TOPOLOGY 

Our new bound m can be turned into a diagnostic 
tool leading to some information on the topology of the 
underlying enzymatic network as follows. By solving the 
transcendental equation m, we find a bound on Mmax 
that depends on the Fano factor F and the affinity A, 



FIG. 5: Minimal effective length of the largest cycle Mmax, 
as given by Eq. m, as a function of the Fano factor F for 
different affinities A. From right to left, the values of the 
affinities are A = 2, A = 5, A = 10 , A = 15, and A —>■ 00 
(which corresponds to the bound Mmax > 1 /F). 


as shown in Fig. O The affinity imposes important con¬ 
straints on the allowed values of F and Mmax- First, 
Eq. (pri) implies F > 2!A. For example, if the free 
energy consumed when one substrate molecule is trans¬ 
formed into a product molecule is 2 A:bT (which implies 
A = 2), the Fano factor F must be larger than 1. Sec¬ 
ond, even far from equilibrium, i.e., for large ratio be¬ 
tween forward and backward transitions rates, the affin¬ 
ity dependent bound on Mmax can represent a consider¬ 
able improvement in relation to Mmax > 1/^- Specif¬ 
ically, if the free energy consumption for transforming 
one molecule is lOfcsT, i.e., A = 10 , a measurement that 
gives a Fano factor F = 0.22 implies that the cycle with 
largest effective length has Mmax > 9, which is almost 
twice the bound obtained for a formally divergent affin¬ 
ity Mmax > 1/F ~ 4.55. We stress that this discussion is 
novel even for the unicyclic case. Furthermore, the size of 
the largest cycle in the network cannot be inferred from 
our bound. The cycle with Mmax does not need to cor¬ 
respond to the cycle with largest N, as demonstrated in 
the second case study from Sec. IIIII 


V. CONGLUSION 

We have found a universal lower bound on the Fano 
factor that takes into account the chemical potential dif¬ 
ference involved in the enzymatic reactions. For unicyclic 
networks, this result, which has been obtained in [l^ . 
is an improvement in relation to the well known bound 
F > n/N, which corresponds to Eq. (1311) with A ^ 00 . 
This limit of divergent affinity happens when some chem¬ 
ical reaction is assumed to be irreversible. Such an as¬ 
sumption typically does not do any harm if one is inter¬ 
ested in the rate of formation of product as a function 
of the substrate concentration, for example. However, in 
principle all chemical reactions are reversible and ther- 
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modynamic consistent models must have a finite affinity. 

Based on several case studies supporting the idea that 
a multycyclic network does not lead to a Fano fac¬ 
tor smaller than the bound associated with the cycle 
with largest effective length M = N/n, we here conjec¬ 
ture that this bound also holds for arbitrary enzymatic 
schemes, with a possibly complex multicyclic network of 
states. Measurements of the Fano factor will now provide 
information on the size of the cycle with largest effective 
length, independent of the details of the underlying en¬ 
zymatic network. Provided the affinity is known, as, e.g., 
in a recent experiment on the Fl-ATPase [l^, our bound 
can represent a considerable improvement in relation to 
Afmax > 1/^j as represented in Fig. [S] 

We have restricted our analysis to enzymatic processes 
for which there is only one chemical potential difference. 
Investigating possible bounds on the Fano factor for the 
case where different substrates and products lead to sev¬ 
eral chemical potential differences would be interesting. 
Finally, on a technical side, we stress that finding an 
algebraic proof of our well-supported conjecture for ar¬ 
bitrary enzymatic networks represents a serious mathe¬ 
matical challenge. 


Appendix A: Expression for the dispersion D 


In this Appendix we explain the method we use to 
calculate the Fano factor, which is due to Koza [l^. We 
consider a general network of states with transition rates 
from i to j denoted by kij. The total number of states 
in the network is Q. We define a Q x Q matrix, which is 
a function of a new variable z, as 


[Mz)y = 


_ J 


k^j exp{zdij) 
— 


iii=j 


(Al) 


The generalized distance dij characterizes how much the 
random variable X changes in a transition from i to j. 
If a transition from i to j involves the enzyme liberating 
a product molecule then dij = 1 and dji = — I. If a 
transition between i and j does not involve a product 
molecule, then dij = 0. As explicit examples we consider 
the two case studies analyzed in section lllll For the model 
in Fig. m L{z) becomes 


-ri 

ya 1 lb 
^-21 w ‘21 

ksi 

0 



12 + 
ki3 

12 “^2 

0 

0 

-r3 

ki2 

k%3 + 

b 

43 

, (A2) 

0 

k24 

Ua _|_ jb 
^34 ^ ^34 

-r4 



b 

12 “ 

yb „-z lb 
'^12^ ’ '■34 

— p ^ 

— /i34e 

]b — 
‘21 ~ 

l.b pZ jb _ 

^21^ ’ M3 ~ 


fc^ge^. For the model in Fig. [31 Liz) is given by 


/ -n 

^21 

k3ie^ 

0 

0 

0 

\ 

ki2 

-r2 

^32 

k42 

fc52e^ 

0 


ki3e~^ 

^23 

-4-3 

0 

^53 

fc63e^ 


0 

^24 

0 

-r4 

^54 

0 


0 

fc25e“^ 

^35 

^45 

-rs 

^65 


V 0 

0 


0 

kse 

-re 

/ 


(A3) 


It can be shown that the maximum eigenvalue of 'L(z), 
which we denote by A(z), is the scaled cumulant generat¬ 
ing function associated with the random variable X [2^ . 
Hence, the following expressions are valid 

V = X' (A4) 

and 

D = X"/2. (A5) 

where the prime denotes derivatives taken at z = 0. The 
explicit calculation of A(z) is usually quite difficult. How¬ 
ever, it is possible to obtain V and D without calculating 
A(z) by using the following method [l^. The character¬ 
istic polynomial related to the matrix L(z) is defined as 

Q 

p{z, y) = det {yl - L(z)) = ^ Cq{z)y'>, (A6) 

9=0 

where Cq{z) are the coefficients of this characteristic 
polynomial. Since A(z) is a root of the polynomial, it 
follows that 

Q 

J2Cq{z)X^z)=0. (A7) 

9=0 

Moreover, for z = 0 the matrix L(z) becomes the stochas¬ 
tic matrix, as given by Eqs. (l22l) and (l28l) for the two 
case studies, which has maximum eigenvalue A(0) = 0. 
Therefore by taking the first derivative with respect to z 
in Eq. (ITtI) we obtain 

^ = (AS) 

and by taking the second derivative we get 

n = (A9) 

where we have used Eqs. (TAdl) and (|A5I) . Hence by 

calculating these coefficients of the characteristic poly¬ 
nomial associated with L(z) we obtain the Fano factor 
F = 2D/V. 


Appendix B: Explicit calculations for a simple 
network 


For illustrative purposes we calculate the Fano factor 
explicitly for the two-state network shown in Fig. jS] We 
assume that both cycles in this network have affinity 


l,a I,b yc yb 


Ua i,b 
^21^12 


Ai2i^l2 


(Bl) 


Assuming that a product is generated in jumps from 1 to 
2 through channels a and c, the modified generator dAH) 
for this model reads 


-ri 

p-2 l^b I l^c 
^12“ ^ ^12 ^ ^12“ 


p 2 _i_ j^b I J^c ^ z 
^2l“ ' ^21 ' ^2l“ 

-T2 


, (B2) 
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a 



FIG. 6: Two-state network used for the explicit calculations. 



FIG. 7: Further examples: four-state network of states. 


where ri = k ^2 + ^12 + ^12 r 2 = k^i + k\i + k^i ■ The 
coefficients of the characteristic polynomial (IA 6 I) associ¬ 
ated with this matrix are ^ 2 ( 2 ) = 1, Ci(z) = ri + r 2 , 
and 


Co(z) = (e-^-I)(e^si-S2), (B3) 

where si = (fc “2 + ^ 12)^21 S 2 = {k^i +^ 21 )^ 12 - Using 
expressions (IA8|) . (IA9I) and F = 2DjJ we obtain the 
following Fano factor, 


Si + S2 2 Si — S2 
Si - S2 (?’i+r2)2’ 


(B4) 


As in the two case studies in the main text the bound 
is saturated when the transition rates are such that 
one cycle with uniform rates dominates. This feature 
can be demonstrated by setting the rates as ^“3 = ^21 = 
exp(yl/ 2 ), ^21 = ^12 = ^12 = fcexp(^/ 2 ), and k^i = k. 

With these rates, the Fano factor (IB4I) becomes 


4e-^/^(l + k) + (e-^ -I-1)(2 + 2k + fc^) 
( e ^ - I )(2 + A :)2 


(B5) 


which is an increasing function oi k > 0. The bound is 
saturated for fc —>• 0 , where 


F = 



(B 6 ) 


Appendix C: Further examples of multicylcic 
networks 

We have also tested our bound expressed in (1311) for 
other multicyclic networks. In this Appendix we discuss 
two examples that we have analyzed extensively. For all 
cases discussed below we verified the bound m with 
two independent checks: a random search in the space of 
transition rates and numerical minimization of the Fano 
factor. 

The first example is the four-state network in Fig. [T] 
This network has three cycles: Cj = (1,2, 3,1), Cjj = 
(1,4,3,1), and C 77 / = (1,2, 3,4,1). The affinities of the 
three-state cycles are „47 = In[/ci 2 fc 23 fc 3 i/(fc 2 ifc 32 fci 3 )] and 


2 



FIG. 8: Further examples: five-state network of states. 


^77 = ln[/(;i 4 /c 43 A: 3 i/(fc 4 ifc 34 fci 3 )]. The affinity of the four- 
state cycle depends on the other two affinities as ^777 = 
Ai — All- The affinity of all cycles must be an integer 
multiplying A = fis — ^J-p, because of the physical con¬ 
straint that there are only one substrate and one prod¬ 
uct. We have tested our bound for the following three 
cases. First, taking Ai = .477 = A and Am = 0 we 
obtain = 3. For this choice the bound is saturated 

if one of the three-state cycles dominate. This case hap¬ 
pens if, for example, we set ki 2 = ^23 = k^i = exp(A./3), 
k 2 i = ^32 = ki 3 = 1, fci 4 = fcexp(A/3), fc 4 i = 1/k, 
ki 3 = exp(^/3)/fc, and ^34 = fc, with A: —>■ 0. Sec¬ 
ond, taking Aj = Am = A and ^77 = 0 we obtain 
Afmax = 4. With this choice the bound is saturated by 
setting the rates as ki 2 = ^23 = ^34 = A: 4 i = exp(A./4), 
^21 = ^32 = ^43 = ^14 = 1, fci 3 = fcexp(^/ 2 ), and 
/cis = k, with A: —>■ 0. Third, we have considered 
the case Ai = —An = A and Am — 2A, for which 
4fmax = 3. The bound is saturated if one of three-state 
cycles dominate, which happens if we choose the rates 
as ki 2 = k23 = ksi = exp(^/3), k 2 i = k32 = ki 3 = 1, 
A:i 4 = A:exp(—2^/3), k^i = l(k, ki 3 = exp(—2^/3)/A;, 
and A :34 = k, with A: —>■ 0. 

The second example is the five-state network in Fig. 
m This network has three cycles: C 7 = (1,2,3,1), 
C 77 = (1,3,4, 5,1), and C 777 = (1, 2,3,4, 5,1). Their 
affinities are Ai = ln[A:i2A:23A:3i/(fc2ifc32^i3)], -477 = 
ln[A:i3A:34A;45A:5i/(A:3iA:43A;54A;i5)], and Am = Ai -f . 477 . 








We have verified our bound for the following two cases. 
First, with Ai = An = A and Am = 2A we obtain 
Afmax = 4. The Fano factor saturates the bound if the 
four-state cycle Cn dominates, e.g., kis = ^34 = ^45 = 
fcsi = exp(^/4), fcai = ^43 = ^54 = fcis = 1, ki 2 = 
fcexp(3^/4), k 2 i = k 23 = 1/fc, and ^32 = k, with fc —> 0. 


Second, with Ai = 8 ^ and Am = —An = AA we obtain 
.^max = 5/4. The bound is then saturated if Cm dom¬ 
inates, e.g., ki 2 = k 23 = ^34 = ^45 = ^51 = exp(4yl/5), 
k 2 i = k32 = *43 = *54 = *15 = 1, *31 = fcexp(32yl/5), 
*13 = k, with * —>■ 0 . 
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